Atomistic Mechanisms of Nonlinear Graphene Growth on Ir Surface 
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As a two-dimensional material, graphene can be naturally obtained via epitaxial growth on a 
suitable substrate. Growth condition optimization usually requires an atomistic level understand- 
ing of the growth mechanism. In this article, we perform a mechanistic study about graphene 
growth on Ir(lll) surface by combining first principles calculations and kinetic Monte Carlo (kMC) 
simulations. Small carbon clusters on the Ir surface are checked first. On terraces, arching chain 
configurations are favorable in energy and they are also of relatively high mobilities. At steps, some 
magic two-dimensional compact structures are identified, which show clear relevance to the nucle- 
ation process. Attachment of carbon species to a graphene edge is then studied. Due to the effect of 
substrate, at some edge sites, atomic carbon attachment becomes thermodynamically unfavorable. 
Graphene growth at these difficult sites has to proceed via cluster attachment, which is the growth 
rate determining step. Based on such an inhomogeneous growth picture, kMC simulations are made 
possible by successfully separating different timescales, and they well reproduce the experimentally 
observed nonlinear kinetics. Different growth rates and nonlinear behaviors are predicted for differ- 
ent graphene orientations, which is consistent with available experimental results. Importantly, as 
a phenomenon originated from lattice mismatch, inhomogeneity revealed in this case is expected to 
be quite universal and it should also make important roles in many other hetero-epitaxial systems. 



I. INTRODUCTION 

Graphene, a monolayer of sp^ hybridized carbon 
atoms, has attracted an intense research interest recently 
due to its unique electronic structure and great applica- 
tion potential^"— Currently, there are several ways avail- 
able to produce graphene. The most elegant method us- 
ing micromechanical exfoliation is ready to produce high 
quality samples,^ but it is difficult to be scaled up. Solu- 
tion based mass production, on the other hand, is chal- 
lenging to reach a high sample quality.— Epitaxial grovifth 
can generate large graphene sample on metal surfaces po- 
tentially with high quality, which makes it the method of 
choice for large scale electronic applications of graphene. 
Various metal surfaces (including Ni, Cu, Ru, Ir, Pd, Co, 
Pt, and etc) have been used to grow graphene<^"— No- 
tice that, on different substrates, the graphene growth 
behavior can be very different.— 

Ir substrate is very attractive for graphene growth, 
where large-scale samples with long-range order and con- 
tinuity can be obtained. Graphene can spread over 
step edges on Ir surface just like a carpet,- and it is 
typically well aligned with the Ir substrate, forming a 
moire patternJ^"— Interestingly, the moire pattern can 
then serve as a template to grow patterned structure of 
metal clusters, with a great potential in nanocatalysis 
and nanomagnetism.-— ^ At the same time, our previous 
study has shown that Ir surface is very rigid compared 
to the soft Cu surfacCfi^iii which makes it a good model 
system to study the growth mechanism of graphene. 

On Ir(lll) surface, several graphene orientations have 
been identified .— Compared to the majority one (RO), 
three other orientations are rotated by 14° (R14), 18.5 ° 
(R18.5), and 30° (R30), respectively. This observation 
indicates that the interaction between graphene and Ir 
substrate is not too weak to be totally orientation inde- 



pendent, and it is also not too strong to only pick up a 
single orientation. A recent angle-resolved photoemission 
spectroscopy (ARPES) experiment has indicated that 
graphene is chemisorbed with its RO orientation while ph- 
ysisorbed with the R30 orientation.^* Therefore, different 
orientations lead to different electronic structures, which 
then provides a flexibility to tune the electronic structure 
of a grown graphene by controlling its orientation. ^8, 

With multiple orientational variants, the growth of 
graphene on Ir surface is complicated. Typically, the RO 
phase nucleates first. Then, at the boundary of a growing 
RO island, one of the three minority orientations (mainly 
R30) can occasionally also nucleate. Once nucleated, the 
R30 phase grows much faster than the RO phase. Most 
interestingly, growth rate of the dominant RO phase has a 
nonlinear (approximately quintic) relationship with the 
concentration of carbon adatoms, which indicates that 
five-membered cluster attachment makes an important 
role in the RO phase growth.— A homogeneous kinetic 
model has been constructed accordingly by fitting exper- 
imental data for a similar surface, Ru(OOOl).^- However, 
the atomistic mechanism of graphene growth on Ir sur- 
face is still unclear. 

In this article, first-principles calculations are carried 
out to study graphene growth on Ir(lll) surface. Adsorp- 
tion, diffusion, and coalescence of small carbon species 
are considered first. Then, nucleation at step sites is 
discussed. Finally, the growth process is studied by at- 
taching small carbon species to a graphene edge described 
with a nanoribbon model. With all these elementary pro- 
cesses understood, the overall picture of graphene growth 
on Ir surface is then obtained, which is characterized by 
inhomogeneous growth and orientation sensitivity. Based 
on this picture, a multiscale kinetic Monte Carlo (kMC) 
model is constructed, which well reproduces the experi- 
mentally observed nonlinear growth behavior. 
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II. COMPUTATIONAL DETAILS 

All calculations were carried out with plane- wave based 
density functional theory (DFT) implementation in the 
Vienna ab initio simulation package (VASP)j^— A 400 
eV kinetic energy cutoff was chosen for plane waves. The 
Perdew-Burke-Ernzerhof (PBE) functional^ was used to 
describe exchange and correlation. K-point sampling 
within the Monhorst-Pack scheme^^ was carefully tested 
to produce well converged results. A four-layer slab 
model was used to describe the Ir(lll) surface, with 
the bottom layer fixed to its optimized bulk geometry 
(the corresponding Ir-Ir bond length is 2.74 A). Repeated 
slabs were separated by more than 10 A to avoid interac- 
tion between each other. Vicinal (322) and (332) surfaces 
were used to model the stepped Ir(lll) surface, which 
contains {100} (A-type) and {111} (B-type) microfacets, 
respectively. Graphene growth front was modeled using 
graphene nanoribbons. For the three minority orienta- 
tions of graphene, only the R30 phase was considered 
since the other two are rarer than R30 by one order of 
magnitude. 

Potential energy of carbon clusters on Ir surface is de- 
fined as 

Ep = {Ec/Ir - Elr)/Nc (1) 

where E(j //^ and Ejr are energies of the absorbed system 
and the clean Ir surface, and Nc is the number of car- 
bon atoms in the cluster. Energy of an isolated carbon 
atom in a big supercell is taken as an energy reference. 
Climbing image nudged elastic band (CI-NEB) method^l 
was used for transition state location and barrier height 
determination. Residual forces were within 0.02 eV/A 
for both geometry optimization and transition state lo- 
cation. 



III. RESULTS AND DISCUSSION 
A. Small carbon species on terrace 

Relative stabilities of different carbon species on 
the surface determines their evolution trend. Here, 
as many as possible configurations of carbon species 
Cn{N=1,. . . ,10), including both chain and compact 
structures, are explored. For a carbon monomer, as 
shown in Figure [1^, the hexagonal close-packed (hep) 
hollow site is the most favorable adsorption site {Ep — 
-7.43 eV), which is 0.26 eV more stable than the face- 
centered cubic (fee) hollow site, agreeing well with previ- 
ous results.?® Carbon atom at a top site is 1.73 eV higher 
in energy than the most stable hep hollow site, with the 
Ir atom under the atop carbon atom being pulled out 
from the surface by about 0.5 A. Bridge-site adsorbed 
carbon atom will be spontaneously relaxed to a hollow 
site. The subsurface octahedral site is also much more 
unstable (1.40 eV) than the hep hollow site, consistent 
with the relatively low carbon solubility in Iri^ 
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FIG. 1. (a) One dimensional chain structures of carbon clus- 
ters on Ir(lll). (b) Some two dimensional compact struc- 
tures. 



For small carbon clusters on the surface, generally the 
most stable configuration is a chain structure, with both 
terminal carbon atoms occupying a hollow site. Longer 
carbon chains prefer to form arches on the surface. Simi- 
lar structures have also been predicted on Ni'^'^ and Cu^ 
surfaces. Another important kind of configurations are 
two dimensional (2D) compact structures (Figure [T|d). 
They are all less stable than one dimensional (ID) chain 
structures, except for Cg. Compact Cq ring is slightly 
more stable (0.06 eV/atom) than its corresponding arch- 
ing chain structure. Although clusters considered in this 
study are relatively small, we can already see a trend 
for compact structures to be dome-like)^ with peripheral 
carbon atoms bonded more strongly to the Ir surface. 

Relative stabilities of carbon species with different sizes 
on Ir(lll) terrace can be compared using the potential 
energy defined in Eq. (1). As a limit case, potential 
energy of graphene on Ir surface is interesting. Its up- 
per bound is estimated by the binding energy of a free 
standing graphene monolayer (-7.99 eV). As shown in 
Figure [U this value is already much lower than those 
of carbon species considered in this study. This is the 
thermodynamic driving force for graphene growth. Car- 
bon monomer is more stable than all small clusters. For 
example, potential energy of carbon dimer is -7.12 eV, 
even higher than that of fcc-hollow carbon. Cg is already 
more stable than fcc-hollow monomer, but still less sta- 
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ble than hcp-hollow carbon. Therefore, concentration of 
small clusters is expected to be much smaller than that 
of carbon monomer on Ir(lll). 




N 



FIG. 2. Potential energy (Ep) of the most stable carbon 
species (Cjv) on Ir(lll) surface and Ir(332) and Ir(322) steps. 
Hollow points indicate 2D compact structures. 




Reaction Coordinate 



Mobility of small carbon clusters on Ir(lll), which is 
mainly determined by the diffusion barrier, is also im- 
portant for graphene growth. The diffusion barrier for 
a carbon monomer from hcp-hollow to fcc-hoUow is 0.71 
eV. As expected, mobility of C2 and C3 is even lower, 
with a diffusion barrier of 0.76 and 0.93 eV, respectively. 
However, for larger carbon clusters with arching chain 
structure, the diffusion barriers are much lower (0.40 eV 
for C5 and 0.54 eV for C4 and Cg). These low barriers are 
a result of the walk-with-legs diffusion mode of the arch- 
ing chain structure (Figure [3]). Although not calculated, 
mobility of two dimensional compact configurations is 
expected to be much lower, since their diffusion usually 
requires to break several bonds simultaneously. When 
mobile carbon species meet, they can coalesce on the 
surface. The energy barrier for the combination of two 
neighboring carbon monomers is 1.44 eV, and it is sim- 
ilar (1.42 eV) for the incorporation of one carbon atom 
into an existing dimer. When adding a carbon atom to 
an arching chain structure of C4, the energy barrier de- 
creases to 0.86 eV. From the energy-barrier point of view, 
formation of small carbon clusters on terrace is generally 
more difficult than their diffusion. 



B. Effects of step edges 

During graphene growth on Ir surface, nucleation 
prefers to occur at step edges 1^ It is thus very desir- 
able to study carbon species adsorption at steps. Both 
A-type and B-type steps are considered by using Ir(322) 
and Ir(332) surfaces, respectively. By subtracting the 
(111) terrace contribution from surface energy, we can 



FIG. 3. Minimum energy path for (a) C5 and (b) Ce car- 
bon chains diffusion on Ir(lll) surface. Inset shows geomet- 
ric configurations before and after each elementary diffusion 
step, which is characterized by hopping of one terminal carbon 
atom between two hollow sites (the walk-with-legs mode). 



compare the stability of these two types of steps. The A- 
type step is more stable than the B-type, with their step 
formation energy as 0.37 and 0.45 eV/A, respectively. 
We study the more stable A-type first. The most stable 
adsorption site of atomic carbon is at the lower edge of 
the step (Figure S]) . The corresponding potential energy 
{Ep=-7.50 eV) is 0.07 eV lower than that of the hcp- 
hollow carbon on terrace, which means that step edge 
is more favorable than terrace. To escape from the step 
edge, a 0.86 eV barrier should be conquered. Dimer and 
trimer also tend to be attracted at the step edge, with a 
potential energy of -7.69 and -7.50 eV, respectively. The 
step edge fully covered with carbon monomers, however, 
is highly unstable (i?p=-7.00 eV). The reason is that car- 
bon atoms will sink into the {100} microfacets, which 
induces a large geometry distortion. The step edge fully 
occupied by dimers has an Ep (-7.36 eV), which is not 
very high and is similar to that of hcp-hollow carbon. 
Therefore, dimer can be trapped at the step edge with a 
relatively high concentration. 

The most stable structure of a carbon trimer at the 
step edge is neither perpendicular nor parallel to the step. 
For larger clusters, a handle-like chain structure perpen- 
dicular to the step is energetically very favorable. In such 
a structure, one terminal atom occupies a bridge-like step 
site just as in the dimer case, while the other terminal 
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FIG. 4. Stable chain structures of carbon species on Ir(322) 
surface. 



atom occupy a hollow site on terrace. Due to the limi- 
tation of the (322) model used for the step, we can not 
study perpendicular carbon chain longer than seven. In 
those cases, we only consider carbon chains along the step 
edge. With significant lattice mismatch between carbon 
and iridium, such chains are typically not well accom- 
modated. Structure of chains with similar length can be 
notably different. This also makes chains along the step 
edge energetically not very favorable. 

Compact structures are generally less stable than the 
chain structures for small carbon clusters at the step 
edge. However, a class of magic compact structures with 
high stability have been identified. The smallest example 
is a C5 cluster (Figure [5]), which has the same energy as 
its corresponding handle-like chain structure {Ep— -7.39 
eV). It can be considered as a structure relaxed from two 
dimers at step edge connected by another carbon atom on 
terrace. Since the step edge can trap many dimers, such 
a structure motif may act as a starting point of graphene 
growth. To check this possibility, we extend this struc- 
ture to two rings (Cg). This structure is already more 
stable than its corresponding handle-like chain structure. 
More importantly, a fragment of zigzag edge with the RO 
orientation appears on terrace in this structure. This 
is consistent with the experimental observations than RO 
phase nucleates at step edges, and it prefers to grow with 
zigzag edgesi^^i^ 



Similarly, this kind of magic structures can be further 
extended to three rings (Cna), which is also the most 
stable structure for clusters with the same size. At the 
same time, instead of extending along the zigzag direc- 
tion, three-ring structure can also be obtained by grow 
an additional hexagon in front of Cg (Cub in Figure [S]). 
This Cub structure can be considered as an initializa- 
tion of graphene growth from a Cg structure, and the 
intermediate between Cg and Cub is the Cio structure. 
Cio, Cna and Cub has almost the same potential energy, 
similar to that of hcp-hollow carbon on terrace. 
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FIG. 5. Stable compact structures of carbon species on 
Ir(322) surface. 

Although it is an elementary unit of graphene, hexagon 
is not stable at the step edge {Ep=-7.08 eV) compared 
with the handle-like chain structure of Cg. It is even 
less stable than hexagon at terrace. Therefore, graphene 
growth is more likely start from the semicircle C5 struc- 
ture instead of forming hexagon at step edges. The most 
stable structure of Cg at step looks like half of the Cio 
structure, which is 0.03 eV per carbon atom more stable 
than the handle-like chain structure. A compact 5-6 ring 
(Cg) becomes more stable than the corresponding chain 
structure along the step edge, but it is still notably less 
stable than neighboring Cg and Ciq. 

Carbon clusters on the Ir(332) surface are similar to 
those on Ir(322) except the missing of magic compact 
structures. Monomer and dimer are very stable at the 
step edges, with a Ep of -7.54 and -7.58 eV, respectively. 
The diffusion barrier of a carbon monomer from a step to 
the lower terrace is 1.05 eV. When atomic carbon fully 
covers the step, its potential energy becomes -7.37 eV. 
Dimer fully covered step gives a potential energy still as 
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low as -7.52 eV. At the B-type step edge, trimer can form 
a straight chain perpendicular to the step, while C4 be- 
comes bent. The most stable structure of Ce and C7 is 
also handle-like chain. However, for C5, the handle-like 
structure has a too large curvature (Figure [6^). There- 
fore, its most stable structure becomes the chain along 
the step edge, which is 0.02 eV lower in potential en- 
ergy. Compact structure becomes more stable than chain 
structure until the cluster size equal to 10, where a fused 
two-hexagon structure is energetically favorable (Ep=- 
7.38 eV). Notice that in this structure there is also a 
fragment of RO oriented zigzag edge. 




FIG. 6. (a) Handle-like structure of C5 perpendicular to the 
Ir(332) step, (b) Chain structure of C5 parallel to the Ir(332) 
step, (c) Compact structure of Cio at the Ir(332) step. 



C. Carbon attachment at graphene edge 

Growth generally proceeds by attaching small carbon 
species on the edges of a growing island. Here, we use 
graphene nanoribbon as a model edge system. To keep 
the supercell reasonably small, C-C bonds are slightly 
stretched (within several percent) to match the Ir lat- 
tice. Consistent with experimentally studied graphene 
islands, which are dome-like with strong edge-substrate 
binding,— nanoribbons also form arched structure in the 
width direction. To study attachment thermodynamics, 
we define an averaged formation energy (Ef) for a spe- 
cific attachment process as 

^ {Ec+R/Ir - En/Jr - Ec'/Ir + Eir)/Nc (2) 

where E(;;_^.ji/ir is the total energy of an absorbed ribbon 
with a carbon species (size Nc) attached, Eft/j^ is the 
energy of the absorbed ribbon before attachment, and 
Ec/ir is the energy of the carbon species absorbed on 
the surface far away before attachment. 

For the RO-phase, we focus on the zigzag edge, since 
it is preferred as observed in scanning tunneling mi- 
croscopy and current imaging tunneling spectroscopy 
experimentsj^i^i Attaching a carbon monomer to a 
zigzag edge (Figure [7]) is energetically unfavorable with 
an Ef equal to 0.58 eV. Since isolated adatom is more 
stable, the attached carbon atom will diffuse away much 
faster than its attachment. Therefore, carbon monomer 
attachment will not contribute to graphene growth at 



such a zigzag edge. In contrary to atomic carbon, four-, 
five-, and six-membered carbon clusters have a driving 
force (about 0.36 eV) to attach to the zigzag edge. With 
endothermic monomer attachment and exothermic clus- 
ter attachment, the experimentally observed nonlinear 
growth seems having a simple thermodynamic mecha- 
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FIG. 7. Optimized structure of a zigzag edge with different 
carbon clusters attached. Ef is marked in eV. 



However, we also note that graphene grown on 
Ir(lll) surface will form a moire structure due to lat- 
tice mismatchfi^ii^ which means that different parts of 
graphene has different registries with the substrate lat- 
tice. For the RO phase, in the atop-type region, all car- 
bon atoms occupy hollow sites. However, in both hcp- 
and fcc-type regions, hollow site carbon is immediately 
neighboring to top site carbon atomsJ^ The model pre- 
sented in Figure [7] thus mainly describes edge fragment 
in the hep- and fcc-type regions, where attached carbon 
atoms occupy top sites. Since top site is energetically 
very unfavorable, attachment of small clusters {N < 4 in 
this case) becomes endothermic in these regions. Notice 
that, for monomer attachment, the attached atom is not 
exactly at a top site in the ribbon model adopted here. 
The structure with an attached atop atom (Figure [HK) 
is a metastable state, which leads to an Ef as high as 
1.98 eV. In atop- type region, the attached carbon atom 
is expected to occupy a hollow site instead of the top 
site. To describe such a behavior, we construct another 
ribbon model. As expected, when the attached carbon 
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atom sits at a hollow site (Figure |8)d), a different ther- 
modynamics is obtained. The atomic carbon attachment 
becomes exothermic with the Ef equal to -0.86 eV. 




FIG. 8. (a) Optimized structure of a carbon monomer at- 
tached at zigzag RO ribbon occupying a (a) top and (b) hol- 
low site, (c) Two carbon monomer and (d) a carbon dimer 
attached R30 zigzag edge on Ir(lll) surface. 

Therefore, attachment of atomic carbon can be en- 
ergetically both favorable or unfavorable depending on 
the growth front configuration. On the kinetic side, the 
energy barrier for carbon monomer attachment for the 
models in Figures [7^ and [Hb are 1.0 and 0.75 eV, re- 
spectively. These barriers are easy to be conquered at 
the experimental temperature. Therefore, as the most 
abundant species on terrace (refer to Figure [5]), car- 
bon monomer's attachment to graphene edge is expected 
to be a process with relatively high possibilities every- 
where. However, for those thermodynamically unfavor- 
able sites, its detachment will be much faster. At those 
sites, graphene growth relies on attachment of much rarer 
carbon clusters. Those parts thus grow much slower, and 
cluster attachment becomes the rate determining step of 
graphene growth. 

To further confirm the above picture, we also perform 
a calculation directly using a large supercell with the 
period in the edge direction matching the exact moire 
unit cell.l- In the direction perpendicular to the edge, 
a narrow ribbon is used, with the other side saturated 
by hydrogen atoms. During geometry optimization, the 
z coordinate of hydrogen atoms is fixed based on the 
experimental mean height of graphene on Ir surface. "^^ 
Along the zigzag edge, eight of the ten sites for carbon 
monomer attachment is energetically favorable. There- 
fore, a structure with eight carbon atoms attached and 
a vacancy at the two top sites (Figure IH^) can be eas- 
ily formed, which is more stable (0.17 eV per attached 
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FIG. 9. (a) graphene ribbon with eight of the ten sites for 
carbon monomer attachment filled and (b) with an additional 
Cs cluster attached. 

carbon atom) than the unattached ribbon plus eight hcp- 
hollow carbon atoms. Although the vacancy is very dif- 
ficult for atomic carbon attachment {Ef=2.11 eV), it is 
well fitted by a C5 cluster (Figure IHb). Such a C5 attach- 
ment process has a.n Ef as large as -0.70 eV. This result 
gives a direct evidence that some parts of the graphene 
edge are difficult for atomic carbon attachment and they 
only grow with carbon clusters supplied. 



D. Kinetic Monte Carlo simulation 

Based on the physical picture from first principles, we 
develop a kMC model to directly reproduce the experi- 
mental nonlinear growth behavior. As already pointed 
out by Zangwill and Vvedensky^SSi a brute-force kMC 
simulation will be impractical for graphene growth on 
Ir(lll). There are mainly two difficulties. First, the 
growth front keeps moving, which requires a huge unit 
cell to get a statistically reliable growth rate. Another 
problem is the disparate rates of different kMC events, 
due to the large density differences for different surface 
species. For example, the gap between densities of car- 
bon monomer and arching chain C5 is about 10^^, which 
stretches the limit of even the most efficient algorithms 
for a direct simulation. To overcome these difficulties, we 
build a multiscale 'growth-front-focused' kMC model. 

Our kMC model is based on a honeycomb lattice. The 
whole surface is divided into four regions: graphene re- 
gion, growth front, diffusion layer, and far field (Figure 
llOb V In the far field, we have an equilibrium of dif- 
ferent A^-sized carbon species. Their density can be ob- 
tained by considering the A^Ci ^ Cat equilibrium"^^ with 
carbon potential energies calculated from first principles. 
Combined with first-principles diffusion barriers, we can 
further get the fiux of difference carbon species across 
the diffusion layer to the graphene growth front. Then, 
we can focus on the carbon species attachment and de- 
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tachment at the growth front, and other processes are 
compacted into effective carbon species fluxes. By doing 
so, we not only reduce the list of kMC events, but also 
extend the grown graphene as long as we need without 
extra cost. 

Rates for the attachment and detachment of different 
carbon species can be calculated, if kinetic data of all 
these events are available from first principles. However, 
due to the effect of the Ir substrate, different sites in the 
honeycomb kMC lattice are not equivalent, which leads 
to a huge number of events. To build a tractable model, 
we make a simplification by considering only two kinds of 
sites. Based on the moire pattern, one kind of sites can 
be easily occupied by carbon monomers, while the other 
is difficult for monomers (see Supporting Information). 
With such a simplification, all required kinetic data is 
calculated or estimated directly from first principles. 









□ kMC J] 




fitting curve ^ 
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0.007 0.010 0.013 0.016 
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FIG. 10. (a) Schematic diagram of the kMC model standing 
on the graphene growth front, (b) Graphene growth rates 
dependence on C monomer concentration obtained by kMC 
simulations. The solid red line is a fitted curve, and the re- 
sulting exponent is 5.17. 

A kMC simulation with the efficient 'growth-front- 
focused' model described above is still prohibited by the 
huge gaps between fluxes of different carbon species. To 
solve this problem, we design a multiscale kMC algorithm 
similar to the nested stochastic simulation algorithm for 
chemical kinetic systems with disparate rates. '^^ In our 
algorithm, the kMC events are grouped by the fluxes of 



the relevant carbon species. Events in different group 
represent the evolvement of the graphene front at differ- 
ent time scales. For example, attachment of C monomers 
represents the fastest time scale. The multiscale kMC 
simulation is run in the following way: 

(1) run a standard 'event-list' kMC algorithm with all 
events included. 

(2) once the graphene front gets stuck, which indi- 
cates that the events on the present time scales don't 
contribute to the graphene growth anymore, turn the 
simulation to the next slower time scale by turning off 
the events in the present time scale group. 

(3) the simulation time scale keeps jumping until the 
graphene grows again. 

(4) once the front of graphene moves, return to step 1. 
With such a specially designed kMC algorithm, we 

run simulations for different C monomer concentrations. 
Convergence of growth rate is carefully tested. The 
obtained relation between growth rate (r) and carbon 
monomer concentration (m) is fitted by r = am'' -I- c. 
As shown in Figure [TUb. we get a h value of 5.17, which 
agree with the experimental value very well. We em- 
phasize that in our kMC simulation, all parameters are 
calculated or estimated from atomic scale, without any 
fitting to macroscale experimental data. Therefore, the 
good agreement between experiment and kMC simula- 
tion suggests that the growth model revealed by our first 
principles calculations has grasped the main physics for 
graphene growth on Ir(lll). 



E. Effects of the Graphene Orientation 

Occasionally, a rotated R30 phase can be nucleated at 
the edge of a growing RO phase. Based on the physical 
picture of the RO phase growth, we can make a predic- 
tion about the growth kinetics of the R30 phase. We start 
from the moire pattern of R30. It is different from that of 
the RO phase, with more homogeneously distributed top 
sitesfi^ In the R30 phase, two top sites are separated by 
at least two hollow or bridge sites. We have shown that 
growth of the RO phase rely on large carbon clusters to 
find their right place where top sites gather. Although 
the moire pattern is not exactly followed at graphene 
edges, the more homogeneous distribution of top sites in 
the R30 phase is expected to remove this stringent large- 
cluster requirement. This will speed up the R30 phase 
growth compared to the RO phase, as exactly observed in 
experiment. It is worth to mention that understanding 
the different growth behaviors for different graphene ori- 
entations will enable us to control the graphene sample 
quality. 

To further study growth of the R30 phase, ribbon mod- 
els are also constructed to describe the growth front. As 
shown in Figure [TTl when a zigzag edge is chosen, carbon 
monomer attachment can be favorable in energy with a 
negative Ef. However, unlike the zigzag edge in the RO 
phase, different edge atoms are not equivalent here. For 
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FIG. 11. Optimized structure of the R30 zigzag edge with 
different carbon species attached on the Ir(lll) surface. 



example, as shown in Figure |SJ3, if the first carbon atom 
occupy an energeticahy favorable hollow site, the next 
one closed to it will occupy a top site, which leads to 
a Ef of 0.83 eV. Similar situation happens for armchair 
edges, carbon monomer attachment at some sites are en- 
ergetically favorable, while it is not preferred at other 
sites. For an energetically favorable monomer attach- 
ment on the zigzag (armchair) edge, the energy barrier 
is 0.99 (1.15) eV. Since it is smaller than that of two 
carbon monomer combination, carbon monomers should 
have the capability of attaching to graphene edges at the 
growth temperature. 

To make a prediction about the kinetics of R30 growth, 
attachment of carbon species with different sizes should 
be studied systematically. As shown in Figure [Til all 
clusters can be attached exothermically at some edge 



sites. For clusters larger than C2, these attachment con- 
figurations already include one or more top sites occu- 
pied. Therefore, if C3 is available, graphene can grow 
over difficult top sites. The possibility of growing over 
top site with even smaller C2 is checked in Figure [8jl, 
which suggests that C2 attachment at top site is still en- 
ergetically favorable, with a negative Ef equal to -0.26 
eV. Since the C2 flux is much larger than those of larger 
clusters, the C2 determined R30 phase growth will be 
much faster than that of RO determined by larger clus- 
ters, and we can predict that the R30 phase growth will 
be more or less quadratic instead of quintic for the RO 
phase. 



IV. CONCLUSION 

Atomistic details of graphene growth on Ir(lll) sur- 
face have been revealed by first principles calculations. 
As building blocks of graphene, small carbon clusters are 
considered first. Arching chain structures are energeti- 
cally favorable on terrace, and they are also of relatively 
high mobilities. Magic compact structures are identi- 
fied at step edges, which are relevant to the nucleation 
process. Graphene nanoribbons are used as models of 
growth front. For the major RO orientation, the growth 
is found to be quite inhomogeneous due to lattice mis- 
match. At some parts of the graphene edge, atomic car- 
bon attachment is thermodynamically prohibited, and 
clusters (size around 5) attachment is required. The 
experimentally observed nonlinear relationship between 
growth rate and carbon monomer concentration is deter- 
mined by such slow cluster attachment steps. A delicate 
kMC model is constructed to directly reproduce the ex- 
perimental growth kinetics. As another result of lattice 
mismatch, graphene growth on Ir(lll) is sensitive to the 
graphene orientation. The rotated R30 phase has less 
inhomogeneity, which leads to a faster growth rate. A 
quadratic growth behavior is predicted for the R30 phase 
growth. Atomistic mechanisms revealed here are essen- 
tial for growth condition optimization in order to control 
graphene orientation and edge morphology, and they may 
also apply to other heteroepitaxial systems with a signif- 
icant lattice mismatch. 

Supporting Information. Potential energy and ge- 
ometry of more carbon species on Ir(lll) and Ir(332), 
more reaction paths of carbon species diffusion and co- 
alescence, cluster attachment for armchair edges, energy 
barriers for carbon atom and clusters attachment, and 
more kMC simulation details. 
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